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ABSTRACT 

Straightforward relationships among Weinstock's 
propagator U A ,the Vlasov propagator U, and the en- 
semble average Vlasov propagator <u> are derived. 
U and <U> are related to the characteristic trajec- 
tories of the Vlasov Equation. 
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PROPAGATORS IN STRONG PLASMA TURBULENCE 

In Dupree's 1 theory of strong plasma turbulence, the fundamental role is played 
by the operator <Jut, t 0 ))>, the average of the Vlasov propagator U(t, t 0 ) over an 
ensemble of plasma realizations. <(U(t,t 0 ))>, as we shall see, can be related to 
various statistical correlations of the turbulent fields. 

Weinstock 2 amplified Dupree's ideas and obtained several formally exact results, 
all involving yet another propagator U A (t, t 0 ). A pair of complicated non-linear 
integro-differential equations implicitly relate U A , U, and <u> in the general 
case . 2 Specific application of the theory 3 is however, limited to the weak 
coupling approximation, where U A is expressed explicitly in terms of (ix) (and 
hence in terms of fluctuation correlations): t j A = <u> to lowest order m a per- 
turbation series in $ F, the amplitude of the fluctuations. 

In this note we first show how in the weak coupling limit U A can be expressed 
straightforwardly in terms of <u> to arbitrary order in b F. Calculation of 
higher order corrections to Dupree's 1 plasma kinetic equation and dispersion 
relation is thus facilitated. We relate U A to <u> in two steps, first relating U A 
to U and then relating U to <U>. Our results are a significant simplification of 
Weinstock's equations. 

We next show the relationship between U (and <U> ) and the characteristic tra- 
jectories of the Vlasov Equation. The value of these latter relationships is that 
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the propagators can be visualized in terms of the Newtonian orbits of Vlasov 
fluid elements and properties of the plasma turbulence. 


Consider an ensemble of Vlasov plasmas. For each realization, c f , the deviation 
of the one particle distribution function from its ensemble average <f> # obeys the 
equation 


frt *■ r-g *(<£”#)•! (1) 


Here < > is the ensemble average and 6 indicates the fluctuation of a quantity, 
from its average value. F is the total force per unit mass on the plasma element 
at the phase space point r , v . 


Weinstock's formal solution to Eq. (1) is 



(2) 


where $f (x,v, t 0 ) is the initial value of S f , L(r) = 6F(r)*B /Bv, and U A (t, t 0 ) 
is defined by Weinstock*s 2 Eq. (8). 


Alternatively we can place the 4j> F • B/ Bv bfy term on the right hand side of 
Eq. (1), recognize B/Bt + y • V + «F> + &£)• B/Bv as the Vlasov operator, 
and iterate with respect to the • BS f/ Bv)> term. The solution obtained in 


this way is 
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U(t,to)SK*Xto) ~ §dVU(t,T) Lfr)<;v(t) 

T tv 
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( 3 ) 


[fdi^, U(r^)LCr^f7(r mH ) - U(r»,t 0 )<W>,s,to)] 

to 


The Vlasov propagator U(t f t 0 ) satisfies the equation 

* if- e u f (<£;<-/£> =o - 1 < 4 > 

The averaging operator A in Eq. (3) averages everything to its right over the 
ensemble of plasmas. If the term <(SF • d8 f /Bv)> were neglected in Eq. (1), 
the solution, Eq. (3), would just be 


io 

After reversing the order of integrations in Eq. (3), changing the variable r n t 
to r, and comparing the resulting form with Eq. (2), we conclude 



< 6 ) 
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Equation 6 is an exact relationship between U A and U. It becomes approximate 
only when the series is truncated. Here A also operates on whatever function 
U A is propagating. 

An equation relating U and <XT> is obtained by averaging Eq. (4) and subtracting 
the resulting equation from Eq. (4) Itself; 


Is *£■ J 7 


<K’M7 (7) 


Equation (7) formally integrates to 


UitM + fdt frt »)) - L(t) (U(T, to)}] (8) 


It is straightforward to iterate Eq. (8) to determine U in terms of <U> to any 
desired order in &F. The term quadratic in U on the right side prohibits us, 
however, from writing U succinctly in terms of <u> to all orders in SF (i.e., 
writing an equation analogous to Eq. 6 between CJ A and U). Iterated through 
Ci(SF) 2 , Eq. 8 is equivalent to Dupree's* Eq. (4.1). 


By combining Eqs. (6) and (8), it is possible to obtain directly a relationship 
between U A and <U> . Through (bF ) 2 , 


jdr, (.'-*) U?,)<0(T h to)') 

fdr, pr A {U it, rj>[(/ - a)/.&)<u c c,, c,j)0 M) <» 

- u z)<m, to))} 


to Z 
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As a check, it is readily shown that with this form of U A inserted into Eq. 2, 
Abf = 0+ S(SF) 3 . This must, of course, be the ease, since * f has by definition 
zero ensemble average* 


We finally indicate the relationship between U and <U> and the characteristic 
trajectories of the Vlasov Equation. If <^F • db f /By^ is neglected in Eq. (1), 
the solution to this equation is 


m,*,*)' fafh'whr]- 


Here x* and v* are solutions to the characteristic equations 




( 11 ) 


(The boundary conditions to be applied are x* (r = t) = x , y* (r = t) = y). 


Next we Taylor expand the x* , v* dependence in Eq. (10) about x, y * Compar- 
ing the result with Eq. (5), we conclude 

U(i, Q -- n? -%]-Vh jV< %) < 12 > 

Operating on an arbitrary function 0(x, v ), U(t, t 0 ) translates the point at which 
0 is evaluated tox (t 0 ), v* (t 0 ), the t = 1 0 phase space coordinates of the plasma 
element located at the point x, v at time t. The trajectory from x* (t Q ), v* (t 0 ) 
to x , v (Eqs. 11) is the exact, fluctuating Vlasov orbit for the element. 



For conciseness of notation we introduce w, the six component phase space 
vector. Thus 






(13) 


In the last form of Eq. (13), we have, following Weinstock 2 , made a cumulant 
expansion 5 . C„ Is the cumulant of fak) •} ] 


We have explicitly written out C 2 in Eq. (13). By integrating the characteristic 
equations, Eqs. (11), A v* (t Q ) can be expressed in terms of the fluctuating fields 
along a particle trajectory. <(U(t, t 0 )^> can thus be represented in terms of 
statistical correlations of the fluctuating field 2 * 6 . One can then further make 
reasonable estimates about the strength of these correlations. 
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